Structure of graphene oxide: thermodynamics versus kinetics 
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Graphene oxide (GO) is an important intermediate to prepare graphene and it is also a versatile material with various 
applications. However, despite its importance, the detailed structure of GO is still unclear. For example, previous 
theoretical studies based on energetics have suggested that hydroxyl chain is an important structural motif of GO, 
which, however, is found to be contrary to nuclear magnetic resonance (NMR) experiment. In this study, we check 
both thermodynamic and kinetic aspects missed previously. First principles thermodynamics gives a free energy based 
stability ordering similar to that based on energetics, and hydroxyl chain is thus thermodynamically still favorable. 
At the same time, by checking the calculated vibrational frequencies, we note that hydroxyl chain structure is also 
inconsistent with infrared experiment. Therefore, kinetics during GO synthesis is expected to make an important role 
in GO structure. Transition state calculations predict large energy barriers between local minima, which suggests 
that experimentally obtained GO has a kinetically constrained structure. 

Recently, an intense research interest has been attracted by graphene oxide (GO).[l| Reduction of GO is a promising 
way to massively product graphene sample. At the same time, many applications of GO have been demonstrated, 
such as electronic device, p-lll chemical catalyst, 0, [To| hydrogen storage, [13, [HI and functional materials. [l3|, EH 
Both to improve the sample quality of graphene from GO reduction and to better utilize GO as a new material, it is 
very desirable to understand the atomic details of GO structure. 

A strong research effort has been devoted to the GO structure study, both theoretically TBI 125^ and experimentally. 
[Bl l26l - l30j Based on nuclear magnetic resonance (NMR) experiment, it is believed that GO mainly has hydroxyl (-OH) 
and epoxy (-0-) groups on the basal plane and carboxylic acid groups (-COOH) at edge sites. [26-281 Hydroxyl and 
epoxy groups are in close proximity, while sp 2 carbon prefers to form small aromatic areas. [22|, [27j More details 
about GO structure, such as possible structural motifs, are still under study. 

Theoretically, many GO structure models based on small supercells have been constructed and asserted by compar- 
ing energies. [TBI, [HI, [TH, Ho|, Hl[ Important GO structure characteristics, such as the proximity of hydroxyl and epoxy 
groups, has been obtained in these studies. However, some controversies between theory and experiment exist. For 
example, although partially oxidized stable GO samples are routinely obtained in experiment, computational ener- 
getics suggests that fully covered GO is the most stable structur e. 1 181 A nother important conclusion from theoretical 
studies is that hydroxyl chain in GO is a very stable structure. [T8l - l21| This result looks reasonable, since hydrogen 
bonds formed in this structure will strongly stabilize the system. However, our previous NMR simulation found that 
hydroxyl chain gives a too large 13 C chemical shift, and it thus should not be an important structural motif of GO. 

A possible reason for these discrepancies is from the thermodynamics: energetically favorable structure may not be 
the thermodynamically most stable one. For example, although GO is synthesized in solution, previous theoretical 
studies did not include solvation effects. Since hydroxyl group can also form hydrogen bond with solvent molecules, the 
stabilization effect of hydroxyl chain may be weakened by taking solution into account. Also, the phonon contribution 
to free energy has not been considered in previous theoretical studies, which may also alter the stability ordering. In 
this study, free energies of different GO structure models have been compared. However, the trend obtained from 
energetics is found to be qualitatively correct in the thermodynamics point of view. Therefore, we also consider the 
kinetics of GO structure evolution. We find that it is difficult to relax the GO structure from one local minimum to 
another. Therefore, GO may exist with a kinetically constrained structure, which explains the controversies between 
theory and experiment well. 

Free energy (G) is calculated with density functional theory (DFT) implemented in the DMOL 3 package [3l|, [32| 
under the generalized gradient approximation [33| 

G = E + AG sol + E ZPE + k B T ^ ln(l - e -^/^ T ) (1) 

i 

where AG so i is the free energy of solvation, and its electrostatic part is calculated using COMSO water solution 
model. [3J, [35j Ezpe is the zero point energy correction, and uji is the phonon frequency at the zone center, which 
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is calculated by finite difference using a large supercell. Minimum-energy pathway for elementary reaction step is 
computed using the nudged elastic band (NEB) method [36] implemented in the Vienna Ab initio Simulation Package 

(vasp). lazllai 

As a benchmark, we first consider isolated epoxy group and hydroxyl group using a 5x5 graphene supercell. With 
a single epoxy group, the optimized structure gives a C-0 bond length of 1.46 A, and the distance between the two 
oxidized C atoms is 1.51 A. The O-H bond length in an isolated hydroxyl group is 0.98 A, and the corresponding C-0 
bond length is 1.49 A. Vibrational frequency of O-H stretch is 3659.1 cm -1 in the isolated hydroxyl group. Geometric 
parameters and vibrational frequency reported here agree well with previous study. [2l[ The calculated AG so i is 
only -0.08 eV for isolated epoxy group, while it is -0.25 eV for isolated hydroxyl group, which is consistent with the 
stronger polarization of hydroxyl group. 




FIG. 1: Different configurations of two epoxy groups. A 5x5 supercell is using in calculation, which is not shown in the figure. 
Carbon is in gray and oxygen is in red. 

To check the trends of energies and free energies ordering for different GO structures, we first consider two epoxy 
groups (20). Four representative configurations including the two with lowest energies [20|, 12JJ are studied (Figured]). 
Their relative energies are listed in Table |TJ la- Id). Structures a and b have almost the same energy, which is 0.10 
and 0.26 eV lower than that of c and d, respectively. This can be understood from a consideration of tension caused 
by epoxy groups. For structures a and b, tensions from both sides will compensate each other to decrease the total 
energy. Aligned epoxy groups prefer to break underlying C-C bonds, pi) as shown in Figure [T]i. The solvation free 
energy varies from -0.08 to -0.18 eV. With the large unit cell used here, zero point energy correction is more than 8 
eV, but their differences are small. Finally, the relative free energies at 300K is similar to relative energies, with the 
same ordering. 

We then study 2 OH models (Figure [2]) with two hydroxyl groups. Their relative energies are listed in Table U (2 a- 
2d). Our calculation shows the 1,2-hydroxyl pair has the lowest energy, which has been widely recognized in previous 
studies. [H, UH Free energies follow the same order as energies, and their relative changes are small. In structures b 
and d, hydrogen bonds are formed. The O-H bonds involving in hydrogen bonds are slightly elongated, leading to two 
softer O-H stretching mode with frequencies about 3450 cm -1 . For other O-H bonds, the corresponding frequencies 
are between 3640 and 3720 cm -1 . 

Considering the proximity of hydroxyl and epoxy groups in real GO structure, we discuss several typical 0+20H 
models (Figure [3]). Their relative energies are listed in Table U (3a-3d). In structure a, two hydroxyl groups are 
attached to carbon atoms directly adjacent to epoxide at the opposite side of the carbon plane, which is the most 
stable. Structure b, with an epoxy group and a neighboring 1,2-hydroxyl pair, has almost the same energy (0.03 eV 
higher). Structures c and d have much higher energies. Calculated free energies have the same trend. Hydrogen bond 
induced softening of O-H stretching mode is also observed similar to the 20H case. 

Besides these isolated-oxidation-group models, we also consider fully oxidized models, which goes to another limit. 
The fully oxidized graphene epoxide (O-full) has two configurations, which are shown in Figure 2^ and b. With a 2x3 
supercell for structure a and a 2x2 supercell for structure b, we have the same number of atoms in both structures. 
The calculated relative energy and relative free energy are still very similar. 

Two models with 100% hydroxyl coverage (OH-full) are constructed. One (Figure |4k) contains hydroxyl chains, 
while the other (Figure 2]i) does not. As shown in Table U (4c-4d) , structure c with hydroxyl chains is 0.13 eV lower 
in energy. Its free energy is also 0.07 eV lower. Strong hydrogen bonds in structure c elongate O-H bonds to about 
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TABLE I: Energy (Eq), solvation free energy (AG so i), zero point correction (Ezpe), and free energy at 300 K ((7(300)) of 
different GO models. Eq and G(300) are listed with the most stable configuration as the reference. All values are in eV. 
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FIG. 2: Different configurations of the two hydroxyl groups. A 5x5 supercell is used in calculation, which is not shown in the 
figure. Carbon is in gray, oxygen is in red, and hydrogen is in white. 



1.0 A and decrease O-H stretching frequencies to 3050-3270 cm -1 . In contrast, O-H stretching in structure d has 
frequencies from 3360 to 3500 cm -1 . 

Structure a in Figure [5] is a stable fully-oxidized GO model (O-OH-full) with hydroxyl chains proposed by Yan et 
al. [HI based on energetics, which is 0.19 eV lower in energy than structure b without hydroxyl chain. When free 
energy is considered instead, structure a is still 0.20 eV more stable than structure b. O-H bonds in hydroxyl chain 
are elongated to about 1.01 A due to the formed hydrogen bonds. Corresponding O-H stretching frequencies are 
generally between 2800-3100 cm -1 . These frequencies for structure b are around 3450 cm -1 . Since the experiment 
value of O-H stretching frequency of GO is about 3400 cm -1 [39|, H(j or higher [4lj, our frequency calculations also 
suggested that hydroxyl chain is not a GO structure motif. 

As another test, partially oxidized structures (O-OH-sp 2 ) with the ratio of C(sp 2 )/C(-0-)/C(-OH) equal to 1:1:1 
are considered (Figure [5b and[5]i). As shown in Table U (5c-5d), structure c with hydroxyl chains is 0.12 eV lower in 
energy than structure d and 0.04 eV lower in free energy. Although the difference becomes smaller when free energy 
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FIG. 3: Different configurations of one epoxy group plus two hydroxyl groups. A 5x5 supercell is used in calculation, which is 
not shown in the figure. Carbon is in gray, oxygen is in red, and hydrogen is in white. 




FIG. 4: (a) and (b) Different configurations of fully oxidized graphene epoxide, (c) and (d) Different configurations of fully 
oxidized GO with hydroxyl groups only. The unit cell used in calculation is marked in red. Carbon is in gray, oxygen is in red, 
and hydrogen is in white. 



is considered, it is safe to conclude that hydroxyl chain is thermodynamically still a very stable structure. The O-H 
stretching frequencies in structure c are among 2800-3100 cm -1 , much smaller than experimental value. Those in 
structure d are around 3580 cm -1 . 

Since hydroxyl chain gives too large chemical shift and too small O-H stretching frequency compared to experimental 
values, while it is still very stable with both solvation and phonon contributions to free energy taken into account, 
it is more likely the formation of hydroxyl chain is prohibited by kinetics during the formation of GO. To check the 
possibility of experimentally obtained GO as a metastable state, we consider the oxidation-group diffusion in GO. 

As shown in Figure [6fci, isolated hydroxyl group is easy to diffuse on pristine graphene surface [24| with an energy 
barrier only 0.27 eV. The C-0 distance is 1.51 A in the initial state and it is elongated to 2.51 A in the transition 
state. When we add an additional neighboring hydroxyl group as shown in Fig [6)3, the diffusion barrier increases to 
1.13 eV. Therefore, the diffusion barrier of hydroxyl group strongly depends on its environment. It is thus important 
to check the hydroxyl mobility in a realistic GO model with a certain epoxy and hydroxyl group concentration. With 
the model shown in Figure [6J3, the diffusion barrier becomes as large as 2.89 eV. This is because that a stable GO 
structure requires a delicate balance between positions of hydroxyl and epoxy groups to minimize stress. Typically, 
moving a hydroxyl group at the same time fixing other groups will lead to a large change in energy. In our case, the 
finial state is 2.57 eV higher in energy than the initial state, which leads to a very high diffusion barrier. 

The diffusion barrier for an isolated epoxy group is 0.74 eV, similar to the value (0.9 eV) obtained by Li et al. 
previously with cluster models. [42| When a neighboring epoxy group is added (Figure [7)3), the diffusion barrier 
becomes 0.40 eV, which also demonstrates a strong environment dependence. Similar to the hydroxyl diffusion case, 



FIG. 5: (a) and (b) Different configurations of fully oxidized GO with both epoxy and hydroxyl groups, (c) and (d) Different 
configurations of partially oxidized GO with both epoxy and hydroxyl groups. The unit cell used in calculation is marked in 
red. Carbon is in gray, oxygen is in red, and hydrogen is in white. 
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FIG. 6: (a) and (b) Energy profile for hydroxyl diffusion. Insets are optimized structure of the initial, transition, and final 
states. A 5x5 supercell is used in calculation, which is not shown in the figure, (c) Initial state of hydroxyl diffusion in a 
partially oxidized GO, and its corresponding (d) final state, where relative energies are marked. Carbon is in gray, oxygen is 
in red, and hydrogen is in white. In (c) and (d), atoms involving in the diffusion in marked in different colors (blue for O and 
green for H). 



in a realistic GO model (Figure Eh), a very large diffusion barrier (2.0 eV) is obtained due to the energy difference 
between the initial and final states. 

Exchange of a hydroxyl group and a neighboring epoxy group can be realized by a H diffusion between them. 
There are two different kinds of transition states, which leads to different final states. The first kind of transition 
states occupy 1,2-sites (Figure [8^) and the second one is on 1,3-sites (Figure [8}}). In the first case, the diffusion 
barrier is very small (0.18 eV), as also reported previously with cluster models. [12| In the second case, we obtain a 
much higher diffusion barrier of 0.88 eV. The longer distance between H and O in the transition state leads to higher 
diffusion barrier. When a realistic GO model is used, the diffusion barrier increase to 1.33 and 1.95 eV for the 1,2- 
and 1,3-cases, respectively. 

According to our calculations, although it is possible to design some delicate structures, in which the barrier of 
oxidation group diffusion is low, generally the diffusion barrier is very high at least in one direction. Therefore, 
a global relaxation of GO structure is very difficult. GO may easily relax to a local minimum, as evidenced also 
by molecular dynamics simulations, [23j but different local minima are expected to be separated by large barriers. 
Experimentally available GO thus typically has a kinetically constrained metastable structure, which explains why 
spectroscopic signal of thermodynamically stable hydroxyl chain structure is not observed in experiment. This result 
can also explain why homogenous phase is formed in partially oxidized GO, although phase separated structure is 
more stable. [20| 

In summary, both thermodynamic and kinetic aspects of GO structure have been considered from first principles. By 
taking the solvation effect and phonon contribution into account, we calculate the free energies of different GO structure 
models. We find that hydroxyl chain is really a very stable structure. However, the calculated vibrational frequencies 
and also previously obtained chemical shifts of this structure is not observed in infrared and NMR experiments. 
Therefore, kinetics during the synthesis of GO must have made an important role in GO structure. We then study 
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FIG. 7: (a) and (b) Energy profile for epoxy group diffusion. Insets are optimized structure of the initial, transition, and final 
states. A 5x5 supercell is used in calculation, which is not shown in the figure, (c) Initial state of epoxy group diffusion in a 
partially oxidized GO, and its corresponding (d) final state, where relative energies are marked. Carbon is in gray, oxygen is 
in red, and hydrogen is in white. In (c) and (d), atoms involving in the diffusion in marked in different colors (blue for O and 
green for H). 
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FIG. 8: (a) and (b) Energy profile for H diffusion from a hydroxy 1 group to a neighboring epoxy group. Insets are optimized 
structure of the initial, transition, and final states. A 5x5 supercell is used in calculation, which is not shown in the figure, 
(c) Initial state of H diffusion with 1,2-transition state in a partially oxidized GO, and its corresponding (d) final state, where 
relative energies are marked, (e) Initial state of H diffusion with 1,3-transition state in a partially oxidized GO, and its 
corresponding (f) final state, where relative energies are marked. Carbon is in gray, oxygen is in red, and hydrogen is in white. 
In (c) and (d), atoms involving in the diffusion in marked in different colors (blue for O and green for H). 



the diffusion of oxidation groups in GO. Typically, there is a very large diffusion barrier. Therefore, experimentally 
available GO will has a kinetically constrained structure. To understand GO structure better, more studies on the 
oxidation process of graphite is highly desirable. 
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